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I. INTRODUCTION 

The unique feature of ultrasonic methods consists in 
the fact that for a solid both an anomalous peak of ul- 
trasonic attenuation and an anomalous change in the ul- 
trasound velocity are observable from experiments^ in 
vicinity of critical point. These phenomena are caused 
by the interaction of low-frequency acoustic oscillations 
with long-lived order parameter fluctuations, which pro- 
duce a random force that perturbs normal acoustic modes 
by means of magnetostrictive spin-phonon interaction. 
In this process, relaxation and fluctuation attenuation 
mechanisms can be distinguished. The relaxation mech- 
anism, which is due to a linear dynamic relationship be- 
tween sound waves and an order parameter^ manifests 
itself only in an ordered phase, where the statistical av- 
erage of the order parameter is nonzero. Since the re- 
laxation of the order parameter near a phase-transition 
point proceeds slowly, this mechanism plays an important 
role in the dissipation of low-frequency acoustic oscilla- 
tions. The fluctuation attenuation mechanism, which is 
determined by a quadratic relation between the deforma- 
tion variables in the Hamiltonian of a system with order 
parameter fluctuations, manifests itself over the entire 
critical temperature range. To date, there are a con- 
siderable number of works that deal with a theoretical 
description of the ultrasonic critical anomalies jSiSiS which 
give an adequate explanation of experimental resultsi^ 

One of the most interesting and important problems 
from both experimental and theoretical viewpoints is the 
study of the influence of structural defects on the dy- 
namical process of sound propagation in solids undergo- 
ing phase transitions. The structural disorder induced by 
impurities or other structural defects plays a key role in 
the behavior of real materials and physical systems. The 
structural disorder breaks the translational symmetry of 
the crystal and thus greatly complicates the theoretical 
description of the material. The influence of disorder is 
particularly important near critical point where behavior 
of a system is characterized by anomalous large response 



on any even weak perturbation. In most investigations 
consideration has been restricted to the case of point- 
like uncorrelated defects^ However, the non-idealities 
of structure cannot be modeled by simple uncorrelated 
defects only. Solids often contain defects of a more com- 
plex structure: linear dislocations, planar grain bound- 
aries, clusters of point-like defects, and so on. 

According to the Harris criterion^ the critical behav- 
ior of only Ising systems is changed by the presence of 
quenched point-like defects. The problem of the influ- 
ence of point-like defects on the critical sound propaga- 
tion in Ising systems has been discussed in paper— with 
the use of a e-expansion in the lowest order of approxi- 
mation. However, o ur p ilot analysis of this phenomenon 
showed that in Ref. Qjl some diagrams which are needed 
for a correct description of the influence of the disorder 
were not considered. Furthermore, our numerous inves- 
tigations of pure and disordered systems performed in 
the two-loop and higher orders of the approximation for 
the three-dimensional (3D) system directly, together with 
the use of methods of series summation, show that the 
predictions made in the lowest order of the approxima- 
tion, especially on the basis of the s expansion, can differ 
strongly from the real critical behavior j 13 i 14 i 15 i 16 i 17 i 18 In 
paper— we have realized a correct field-theoretical de- 
scription of dynamical effects of the influence of point- 
like defects on acoustic anomalies near temperature of 
the second-order phase transition for 3D Ising systems in 
the two-loop approximation with consideration of only 
the fluctuation attenuation mechanism. In paper o 20 i 21 
we have performed, the field-theoretical description of the 
influence of point-like defects on critical acoustic anoma- 
lies with allowance for both fluctuation and relaxation 
attenuation mechanisms. 

The temperature and frequency dependencies of the 
scaling functions for ultrasonic attenuation and sound 
velocity dispersion were calculated in paper— the two- 
loop approximation with the use of the resummation of 
asymptotic series technique for pure and disordered sys- 
tems, and their asymptotic behavior in hydrodynamic 
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and critical regions were distinguished. It was shown 
that the presence of structural defects causes a stronger 
temperature and frequency dependencies of the acous- 
tic characteristics in the critical region and increase in 
the critical anomalous peak of ultrasonic attenuation in 
comparison with pure system. 

However, the question about influence of long-range 
(LR) correlation effects realizing in defects with complex 
structure on critical behavior of systems was investigated 
noticeably weaker, and this question has not been stud- 
ied entirely in respect to display of correlation effects in 
the critical anomalous properties of sound propagation 
in solids. Different models of structural disorder have 
arisen as an attempt to describe complicated defects. In 
this paper, we concentrate on the model of Weinrib and 
Halperin (WH)S with the so-called LR-correlated dis- 
order when pair correlation function for point-like de- 
fects g(x — y) falls off with distance as a power law 
g(x — y) ~ \x — y\~ a ■ Weinrib and Halperin showed that 
for a > d LR correlations are irrelevant and the usual 
short-range (SR) Harris criterionii 2 — dv$ = ctQ > 
of the effect of point-like uncorrelated defects is realized, 
where d is the spatial dimension, and i>a and ao are the 
correlation-length and the specific-heat exponents of the 
pure system. For a < d the extended criterion 2—avo > 
of the effect of disorder on the critical behavior was estab- 
lished. As a result, a wider class of disordered systems, 
not only the 3D diluted Ising model with point-like uncor- 
related defects, can be characterized by a different type 
of critical behavior. So, for a < d a LR-disorder sta- 
ble fixed point (FP) of the renormalization group (RG) 
recursion relations for systems with a number of compo- 
nents of the order parameter m > 2 was discovered. The 
critical exponents were calculated in the one-loop approx- 
imation using a double expansion in e = 4 — d <C I and 
6 = 4 — a -c f ■ The correlation- length exponent was eval- 
uated in this linear approximation as v = 2/ a and it was 
argued that this scaling relation is exact and also holds 
in higher order approximation. In the case m = 1 the ac- 
cidental degeneracy of the recursion relations in the one- 
loop approximation did not permit to find LR-disorder 
stable FP. Korzhenevskii et al^ proved the existence of 
the LR-disorder stable FP for the one-component WH 
model and also found characteristics of this type of crit- 
ical behavior. 

The results for WH model with LR correlated de- 
fects received based on the e, 5 expansio n 22 i 23 ' 24 i 25 ' 26 were 
questioned in our paper^ where a renormalization anal- 
ysis of scaling functions was carried out directly for the 
3D systems in the two-loop approximation with the val- 
ues of a in the range 2 < a < 3, and the FP's correspond- 
ing to stability of various types of critical behavior were 
identified. The static and dynamic critical exponents in 
the two-loop approximation were calculated with the use 
of the Pade-Borel summation technique. The results ob- 
tained in Ref. [l6| essentially differ from the results eval- 
uated by a double e, S - expansion. The comparison of 
calculated the exponent v values and ratio 2/ a showed 



the violation of the relation v = 2/ a, assumed in Ref. [22] 
as exact. Monte-Carlo simulations of the short-time dy- 
namic behavior were carried out in paper— for 3D Ising 
and XY models with LR correlated disorder at critical- 
ity, in the case corresponding to linear defects with a = 2. 
It was shown that the obtained values of the static and 
dynamic critical exponents are in a good agreement with 
our results of the field-theoretic description of the critical 
behavior of these models in Ref. 

The models with LR-correlated quenched defects have 
both theoretical interest due to the possibility of pre- 
dicting types of critical behavior in disordered systems 
and experimental interest due to the possibility of real- 
izing LR-correlated defects in the orientational glasses, 28 
polymers^ and disordered solids containing fractal-like 
defects^, dislocations near the sample surface^ or 4 He 
in porous media such as silica aerogels and xerogels^L For 
this purpose, in this work we performed a field-theoretical 
description of the effect of LR-correlated quenched de- 
fects on the anomalous critical ultrasound attenuation in 
3D systems with allowance for both the fluctuation and 
relaxation attenuation mechanisms without using the s, 8 
expansion. We considered in this paper only Ising-like 
systems that gives to compare the effect of LR-correlated 
defects with point-like uncorrelated defects on the critical 
ultrasound attenuation. 

In Sec. II, we describe the model of compressible Ising- 
like systems with LR-correlated disorder and set of dy- 
namical equations which are used for description of the 
critical sound propagation in solids. We also give in 
Sec. II the results of RG expansions in a two-loop ap- 
proximation for ultrasound scaling functions, and (3 and 
7 functions of the Callan-Symanzik RG equations. The 
FPs corresponding to the stability of various types of crit- 
ical behavior and the critical exponents are determined 
in Sec. III. Section IV is devoted to analysis of the ultra- 
sound critical characteristics and discussion of the main 
results. 



II. MODEL AND DESCRIPTION OF 
DYNAMICAL EQUATIONS. RG EXPANSIONS 

For description of the critical sound propagation in 
solids the model of compressible media is used. Inter- 
action of the order parameter with elastic deformations 
plays a significant role in the critical behavior of the com- 
pressible system. It was shown for the first time in Ref. [3^ 
that the critical behavior of a system with elastic degrees 
of freedom is unstable with respect to the connection of 
the order parameter with acoustic modes and a first-order 
phase transition is realized. However, the conclusions of 
Ref. l32 j ar e only valid at low pressures. It was shown 
in Ref. 33 that in the range of high pressures, beginning 
from a threshold value of pressure, the deformational ef- 
fects induced by the external pressure lead to a change 
in type of the phase transition. 

The Hamiltonian of a disordered compressible Ising 
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model can be written as 

H = H c i + H p + H ln t + H mp , 



(1) 



consisting of four contributions. 

The contribution H e \ in ([1]) of the deformation degrees 
of freedom is determined as 
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where u a p{x) are components of the strain tensor and 
C^- are the elastic moduli. The use of an isotropy ap- 
proximation for H e i is caused by the fact that in the 
critical region the behavior of Hamiltonian parameters 
of anisotropic system is determined by an isotropic FP 
of RG transformations, while the anisotropy effects are 
negligible*^ This phenomenon of the increase in system 
symmetry in the critical point is well known as asymp- 
totic symmetry. 

H op is a magnetic part in the appropriate 
GinzburgLandau- Wilson form: 



d d x 



IroS' + l^sf + luoS^ 



(3) 



where S(x) is the Ising-field variable which is associated 
with the spin order parameter, uq is a positive constant 
and To ~ (T — To c ) /Tq c with the mean-field phase tran- 
sition temperature T 0c . 

The term H m t describes the spin-elastic interaction 
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which is bilinear in the spin order parameter and linear 
in deformations, go is the bare coupling constant. 

The term -f/i mp of the Hamiltonian determines the in- 
fluence of disorder and it is considered in the following 
form: 



H imp = Jd d x [At(x)S 2 ]+ Jd d 



, (5) 



where the random Gaussian variables At (a;) and h(x) 
determine the local transition temperature fluctuations 
and induced random stress, respectively. The second mo- 
ment of distribution of the At(x) fluctuations g(x — y) = 
((AT(x)AT(y))) ~ \x — y\~ a characterizes the LR corre- 
lation effects realizing in defects with complex structure. 
The second moment of distribution of the h(x) fluctua- 
tions C(x — y) = ((h(x)h(y))) presets a renormalization 
of the elastic moduli and the coupling constant in the 
spin-elastic interaction. 



To perform calculations, it is convenient to use the 
Fourier components of the deformation variables in the 
form: 

Uaf} = U^l + V~ 1/2 V] U ap(l) ex P (iqx) 7 (6) 



where q is the wavevector, V is the volume, u^l is an 
uniform part of the deformation tensor, and u a p(q) = 
i/2 [q a up + qpu a ]- Then the normal-mode expansion is 
introduced as 

u 0?) = e A(e)Qg,A, 

A 

with the normal coordinate Q q ,\ and polarization vector 
e\(q). We carry out the integration in the partition func- 
tion with respect to the nondiagonal components of the 
uniform part of the deformation tensor ufl, which are 
insignificant for the critical behavior of the system in an 
elastically isotropic medium. After all of the transforma- 
tions, the effective Hamiltonian of the system has become 
in the form of a functional for the spin order parameter 
S (q) and the normal coordinates Q\ (q): 

H=\J d d 9 (ro + q 2 ) S q S- q + J d d qqh q Q- q 
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The Fourier transformation of the correlation function for 
distribution defects in system g(x) ~ |x| _a gives g(k) = 
vo + wok a ~ d for small k. g(k) must be positive definite, 
therefore, if a > d, then the wo term is irrelevant, vo > 0, 
and the Hamiltonian (7) describes the model with SR- 
correlated defects, while if a < d, then the wq term is 
dominant at small k and wq > and characterizes a 
large influence of LR-correlated defects on the critical 
behavior. 

The critical dynamics of the compressible system in 
the relaxational regime can be described by the Langevin 
equations^ for the spin order parameter S q and phonon 
normal coordinates Q q .y. 



S q — —T 



Q q ,\ 
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where Tq and Dq are the bare kinetic coefficients, £ g (x, t) 
and rj q (x,t) are Gaussian white noises, and h$ and hq 
are the fields thermodynamically conjugated to the spin 
and deformation variables, respectively. 

The quantities of interest are the response functions 
G(q,u) and D{q 1 oj) of the spin and deformation vari- 
ables, respectively. It can be obtained by linearization in 
correspondent fields that 



D(q,u>) 
G(q,u>) 



*[(Q q ,u>,\)] 



Sh Q 



[(Q q , u ,xQ- q ,- u ,x)] , (9) 
= [{S 9lU S- 9 ,- u )] , (10) 



where (...) denotes averaging over Gaussian white noises, 
[...] denotes averaging over random fields Ar_ g and h q 
that are specified by structural defects. 

The response functions may be expressed in terms of 
self-energy parts: 

G-\q,u) = Go\q,uj) + U(q,uj), (11) 
D-\q,w) = D- 1 (q,Cj) + j:(q,Lj), 

where the free response functions Go (q, u>) and Dq (q, u>) 
have the forms 

Do(q, ui, A) = 1 /(w 2 — aq 2 — iujD^q 2 } , 
G (q,Lo) = l/[icj/r +(r Q +q 2 )] . 

In the low-temperature phase, the response functions 
contain an additional relaxation contribution caused by 
the presence in the spin density 

S q = MS q>0 + ifiq (12) 

of an additional part connected with the magnetization 



M 



0, 



B\T~T C \ P , T <T C 



(13) 



where B is the phenomenological relaxation parameter 
and ip q is the fluctuation part of the order parameter. 

The characteristics of the critical sound propagation 
are defined by means of the response function 3 ^ D(q,ui). 
Thus, the coefficient of ultrasonic attenuation is deter- 
mined through the imaginary part of H(q,u): 



a(u>, t) ~ u> lm£(0, u>), 



(14) 



and the sound velocity dispersion is expressed through 
its real part, 



c 2 (lu, t) - c 2 (0, t) - Re [E(0, u>) - £(0, 0)] 



(15) 



We calculated E(g,w) in the two-loop approximation. 
The diagrammatic representation of T,(g, u>) is presented 
in Fig. [Q 

The Feynman diagrams involve momentum integra- 
tions in dimension d (in our case, d = 3). Near the 
critical point the correlation length £ increases infinitely. 



When <C A, where A is a cut-off in the momentum- 
space integrals (the cut-off A serves to specify the basic 
length scale) , the vertex functions are expected to display 
an asymptotic scaling behavior for wavenumbers q <C A. 
Therefore, one is lead to consider the vertex functions in 
the limit A — + oo . The use of the RG scheme removes all 
divergences which arise in thermodynamic variables and 
kinetic coefficients in this limit. 

To calculate the attenuation coefficient ([13)1 and the 
ultrasound velocity dispersion (JTSJ) and to eliminate the 
divergences in Y,(q,u>) at q — * 0, we used the matching 
method which was introduced for statics in Ref. [3^ and 
then generalized for critical dynamics in Ref. I3SL First, 
we use the dynamical scaling property of the response 
function 

D (q, co, r) - e^ l D (qe\ (uj/T ) e zl ,re^) , (16) 

and then we calculate the right-hand side of this equation 
for some value I* = I, where the arguments do not all 
vanish simultaneously^ The choice of I* is determined 
by the condition 



(w/r ) e 2 
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= 1, 



(17) 



which permits to achieve the main purpose of the RG 
transformations, namely, to find a relation between the 
behavior of the system in the precritical regime at a low 
value of the reduced temperature r and the behavior of 
the system in a regime far from the critical mode, i.e., 
without divergences in E(q,w). As it was demonstrated 
in Ref. [H, the matching condition (| 1 T[) provides an in- 
frared cut-off for all diverging quantities. The particular 
form of the matching condition (fT7|) permits an explicit 
solution for I* in the form of a functional dependence 
on lj and r, which is specified by the dynamic critical 
exponent z and the static critical exponent v for the cor- 
relation length, 



1 + (y/2) 



-1/4 

4/2 =r-"F(y), 



(18) 



where the abbreviation y — lot~ zv /To is introduced as 
the argument of the F(y) function. 

The response function D [qe l , (uj/Tq) e zl , re l ^ u ] on the 
right-hand side of (fT6|) is represented by Dyson Fq. (fTTj) 
and the following scaling relationships for the imaginary 
and real components of self-energy part are valid: 36 



ImE(tj) 



lmS,{uje zl ) 



(19) 



Re[E(0,w) - £(0,0)] = 

= e Za / !y Re[E(0,we z/ )-E(0,0)] , (20) 

where a is the exponent for heat capacity. It may be 
argued 3 ^ that condition ([T7]) with the well-known expres- 
sion for the susceptibility!? provides an infrared cut-off 
for all divergent values. 
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FIG. 1: Diagrammatic representation of E(g,o;) in a two-loop approximation, (a) Solid line corresponds to Go(q, w); (b) lines 
with a cross, to Co(q, to) — 2F _1 ^(cj/ro) 2 + (q z + t ) 2 J ; (c) vertex with dashed line, to v = [(At) 2 ] ; (d) vertex with waved 
line, to w, which characterized long-range correlation of disorder with wave vector q a ~ d \ and (e) line, to relaxation insertion 
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It was shown in later theoretical works^S that in 
asymptotic regions with r — > and lo — ► the coef- 
ficient of attenuation and the sound velocity dispersion 
are described by a simple scaling functions of the vari- 
able y only. The experimental investigations performed 
on the crystals Gd 8 and MnP 9 confirmed the validity of 
the concepts of dynamical scaling. 

In the asymptotic limit (t — > 0, u) — > 0), the expression 
for the imaginary part of £(w,t) can be defined by a 
scaling function <j){y) 



ImS( W )/ W =r- a -^(|/), 



(21) 



and the expression for the real part of E(w,r) can be 
determined using another scaling function f(y) 

Re [E(0, w) - E(0, 0)] = r~ a [f(y) - /(0)] . (22) 

The substitution of e l from Eq. (fT5)) into the right-hand 



side of the expressions (jl9|) and ((20)) allows to calculate 
the 4>(y) and /(y) scaling functions. 

The dynamic scaling functions calculated in the two- 
loop approximation have the forms: 



g* 2 T Q F a / iy + 1 / 2l/ - 
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where g*, u*, v* and w* are values of the interaction 
vertices in the FP of the RG transformations that cor- 
responds to the critical behavior of the disordered com- 
pressible Ising model,~ generalized in this paper to the 
case with LR-correlated disorder. The terms in Eqs. (|2"3"|) 
and (|24[) that are proportional to M 2 describe the relax- 
ation contribution for the scaling functions of the atten- 
uation coefficient and the sound velocity dispersion. 

In accordance with our paper ^ the values of g* , u* , 
v*, and w* can be obtained from the RG considerations of 
the replicated Hamiltonian for the disordered compress- 
ible Ising model in which we averaged over deformation 
variables. The final form of this Hamiltonian can be pre- 
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TABLE I: Coefficients for the /3- functions in Eq. [26] 



a 


bi 


b 2 


bz 


b 4 


65 


be 


67 


3.00 


1.851852 


9.703 704 


3.425 924 


6.851852 


1.000 000 


2.000 000 


2.666 664 


2.90 


1.751381 


9.149 428 


3.047 432 


6.456 936 


1.632 512 


1.836 920 


2.376 632 


2.80 


1.662 830 


8.671819 


2.747 992 


6.117964 


2.148 800 


1.696 968 


2.190 416 


2.70 


1.584 520 


8.260 292 


2.508 396 


5.827 228 


2.584 256 


1.576 244 


2.086 520 


2.60 


1.515 077 


7.906 550 


2.313 792 


5.578 740 


2.967 744 


1.471 748 


2.053 896 


2.50 


1.453 357 


7.604 029 


2.151672 


5.367 788 


3.325 888 


1.381 120 


2.089 720 


2.40 


1.398 383 


7.347 527 


2.010 060 


5.190 684 


4.022 776 


1.302 500 


2.199 488 


2.30 


1.349 314 


7.132 943 


1.875 288 


5.044 576 


4.904458 


1.234 384 


2.396 896 


2.20 


1.305 402 


6.957111 


1.728 540 


4.927324 


5.918 144 


1.175 580 


2.708 976 


2.10 


1.265 968 


6.817 670 


1.539 252 


4.837412 


7.035 808 


1.125 124 


3.183 336 


2.00 


1.230 378 


6.713 001 


1.250 616 


4.773 916 


8.216 616 


1.082 260 


3.905 832 
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where g(qi + q%) — vo + wo|qi + q2\ a ~ d is the Fourier 
transformation of the correlation function for the 
quenched random distribution of defects in system. The 
Hamiltonian in (|25[) is a functional of n replications of 
the original order parameter. The properties of the orig- 
inal disordered system after the RG transformations of 
Hr are obtained in the replica number limit n — > 0. 

The effective parameter of interaction uq = uq — 2g^ 
arising in the Hamiltonian (|25| from influence of spin- 
elastic interaction determined by the parameter go can 
take both positive and negative values. As a result, the 
Hamiltonian in (j2l)|) describes a phase transitions both 
second order and first order. The tricritical behavior is 
realized in the system for Uq = 0. Furthermore, as it 
was shown in Ref. 10. the dependence of the effective 
interaction in (|2"5|) determined by the vertex g^ — p, (P) 
on outside pressure P can lead to the phase transition of 
second order for P > P t and lead to the phase transition 
of first order for P < P t , where Pt is a tricritical value of 
outside pressure. But we restrict ourselves in this paper 
by consideration of the critical behavior only. 

As is known, in the field-theory approach^ the asymp- 
totic critical behavior of systems in the fluctuation region 
is determined by the Callan-Symanzik RG equation for 
the vertex parts of the irreducible Greens functions. To 
calculate the (3 functions and the critical exponents as 
scaling 7 functions of the renormalized interaction ver- 
' — ft, v, and w appearing in the RG equa- 



tices u, g 2 

tion, we used the method based on the Feynman diagram 
technique and the renormalization procedure^ 



As a result, we obtained the (3 and 7 functions in the 
two-loop approximation in the form of the expansion se- 
ries in renormalized vertices tt, <? 2 , v, and w. We list here 
the resulting expansions: 

flu(u, v, w) = — u + 9u 2 — 2Auv — 16(3/i — f2)ww 
308 o 5920. 2 1664„ 2 

" — u ~^m uv + — uv 

+ 16(&i + b2)u 2 w — QAb^iiw 2 — QAbiUvw, 
Pg(u, g, v,w)= - g 2 + 2<7 4 + 6ug 2 - 8vg 2 - Wfiwg 2 

92 _o_o 224. _ 2 736 2 _ 2 

— — u q H uvq V 5 

9 y 9 y 27 y 

736, 2 . 2 112, . . 2 736, _ 2 

— -7fjr b 5 w 9 + -^-b 6 uwg - —b 6 vwg , 

(3 v (u, v, w) = - v - 16v 2 - I6/3W 2 - 16(3/i - f 2 )vw 

3040 o „„, , 800. n 

+ 6uv v d - 2>2b 7 w A -I uv 

27 9 

92 

+ 12b 8 uw 2 - 32b 9 vw 2 - yjj 2 t) 

— 326iof 2 w + 2ibnuvw, 
w (u,v,w) = -(4 - a)w - 16(/i - f2)w 2 - 8vw 

+ 6uw + 32&12W 3 — 32bi3vw 2 — 27~ v2w 

92. 2 „, . 2 368 . 

— — u w + I2614UW H uvw, 

9 9 



^{u,v,w)= 8f 2 w + ^u 2 + ^v 2 + 8aw 2 
32 

— ~~^~UV — 6C2UW + 8C2VW, 

7^2({(, g, v, w) = — 3u + Av — 2g 2 + 8f\w + 6u 2 + 16w 2 



+ 2g + 8C3W + 8C4VW — 6C4UW — 2Auv, 



h = 
h = 
h = 



(a-2)(a-4) 

2sin(7ra/2) ' 
(a-2)(o-3)(a-4) 

487rsin[7r(a/2 - 1)] ' 
(2a- 5) (2a -7) 



2sin[7r(a-3/2)]' 
where the coefficients bi and c, for different values of pa- 



(26) 
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TABLE II: Coefficients for the /^-functions in Eq. [26] (the continuation of Table HJ. 



a 


bs 


b s 


bio 


fen 


fel2 


fel3 


fel4 


3.00 


4.000 000 


8.851 848 


9.703 704 


5.703 704 


-0.851848 


1.703 704 


3.407 408 


2.90 


3.847 968 


8.018 848 


9.149 424 


5.384 984 


-0.630 048 


1.471 544 


2.943 084 


2.80 


3.777 176 


7.384 568 


8.671816 


5.107 714 


-0.442 952 


1.259 744 


2.519 492 


2.70 


3 785 160 


6 909 368 


8 260 288 


4.866 124 


—0 278 072 


1 063 416 


2.126 828 


2.60 


3.875 560 


6.565 360 


7.906 552 


4.655 620 


-0.125 208 


0.878 440 


1.756 888 


2.50 


4.059 204 


6.332 944 


7.604 032 


4.472 486 


0.024 984 


0.701 240 


1.402 484 


2.40 


4.356 584 


6.198 408 


7.347528 


4.313 686 


0.182 264 


0.528 528 


1.057056 


2.30 


4.802 876 


6.152 008 


7.132 944 


4.176 734 


0.358 480 


0.357160 


0.714 324 


2.20 


5.457 744 


6.185 952 


6.957112 


4.059 572 


0.570112 


0.184 008 


0.368 020 


2.10 


6.425 424 


6.291216 


6.817672 


3.960 514 


0.842 960 


0.005 800 


0.011596 


2.00 


7.899 532 


6.451 000 


6.713 000 


3.878 172 


1.222 400 


-0.181032 


-0.362 064 



TABLE III: Coefficients for the 7-functions in Eqs. ([26)1 . 



a 


Cl 


C2 


C3 


C4 


3.00 


0.296 288 


0.592 592 


2.080 000 


4.000 000 


2.90 


0.423 680 


0.730 760 


1.683 744 


3.673 840 


2.80 


0.539 712 


0.874 448 


1.425 600 


3.393 936 


2.70 


0.651712 


1.025 656 


1.207 872 


3.152 488 


2.60 


0.765 696 


1.186 608 


1.016128 


2.943 496 


2.50 


0.887 040 


1.359 760 


0.837056 


2.762 240 


2.40 


1.021088 


1.547944 


0.656 544 


2.605 000 


2.30 


1.173 728 


1.754 448 


0.456 736 


2.468 768 


2.20 


1.351808 


1.983 136 


0.211584 


2.351 160 


2.10 


1.563 808 


2.238 656 


-0.122112 


2.250 248 


2.00 


1.820 576 


2.526 584 


-0.624 224 


2.164 520 



rameter a in the range 2 < a < 3 are given in Tables ITl llIII 
In the series (|26l) it was used a standard change in vari- 
ables u — > u/J, g 2 — * g 2 / J , v — > v/J, and w — ► w/J, 
where J = J d d q/(q 2 + l) 2 is the one-loop integral. 

We should like note that WH model with a = 3 cor- 
responds to a system with point-like defects. The case 
with a = 2 corresponds to a system of straight lines of 
impurities or straight dislocation lines of random orien- 
tations in a sample. The cases with noninteger a with 
2 < a < 3 are treated in terms of a fractal dimension 
of impurities. We think that the features of the critical 
behavior of disordered systems with 2 < a < 3 can be 
displayed in real highly disordered systems when concen- 
tration of defects is sufficiently high that effective long- 
range correlations can occur because of elastic interaction 
of defects. For values of correlation parameter a less than 
1.5, the unavoidable divergencies appear in diagrams for 
RG j3- and 7-functions in framework of theoretical-field 
approach with fixed dimension d — 3, which is used in 
this paper. 

III. FP'S AND VARIOUS TYPES OF CRITICAL 
BEHAVIOR. CRITICAL EXPONENTS 

The nature of the critical behavior is determined by the 
existence of a stable FP satisfying the system of equations 

0i(u*,g*,v*,w*)=O (i = 1,2,3,4). (27) 



It is well known that perturbation series are asymptotic 
series, and that the vertices describing the interaction of 
the order parameter fluctuations in the fluctuating region 
t — > are large enough so that expressions (|26p cannot 
be used directly. For this reason, to extract the required 
physical information from the obtained expressions, we 
employed the Pade-Borel approximation of the summa- 
tion of asymptotic series extended to the multiparameter 
caseJ£ The direct and the inverse Borel transformations 
for the multiparameter case have the form 



f(ui,U 2 , U 3 ,U4) = 2J c i]lkU\ 
i,j,l,k 



l u J 2 u l 3 u k 4 = 



e 1 B(uit,ii2t,U3t,U4t)dt, (28) 



S, \ \ ^ Cijlk ill k 

{Ui,U 2 ,U 3 ,U4) = > 7— . , U\U J 2 U a U%. 



A series in the auxiliary variable 9 is introduced for an- 
alytical continuation of the Borel transform of the func- 
tion: 



B(ui,u 2 ,u 3 ,u 4 ,8) = (29) 

00 k k—i k—i—j 

fe=0 i=0 1=0 1=0 



Cij. I, k—i—j— I i j I k — i — j — I 

^ — — fc , u\u J 2 u 3 u 4 
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TABLE IV: Stable fixed points and critical exponents for the compressible Ising model with LR-correlated disorder. 



a 


u* 


r 2 


V* 


v* + w* 


V 


Q 


P 


z 


3.01 


0.26482 


0.02043 


0.03448 


0.03448 


0.6896 


— 0.0687 


0.3561 


2.1712 


3.00 


0.26482 


0.02043 


0.01393 


0.03448 


0.6896 


— 0.0687 


0.3561 


2.1712 


2.90 


0.28867 


0.06021 


0.01993 


0.04257 


0.7387 


-0.2160 


0.3806 


2.2120 


2.80 


0.30881 


0.05443 


0.02510 


0.04894 


0.7411 


-0.2232 


0.3805 


2.2486 


2.70 


0.32670 


0.04708 


0.02968 


0.05422 


0.7402 


-0.2205 


0.3785 


2.2837 


2.60 


0.34294 


0.03847 


0.03380 


0.05870 


0.7369 


—0.2106 


0.3749 


2.3184 


2.50 


0.35776 


0.02899 


0.03752 


0.06248 


0.7320 


-0.1958 


0.3703 


2.3532 


2.40 


0.37120 


0.02281 


0.04086 


0.06562 


0.7297 


-0.1890 


0.3668 


2.3879 


2.30 


0.38313 


0.01772 


0.04380 


0.06811 


0.7278 


-0.1833 


0.3634 


2.4215 


2.20 


0.39322 


0.01281 


0.04631 


0.06989 


0.7254 


-0.1761 


0.3597 


2.4524 


2.10 


0.40090 


0.00762 


0.04829 


0.07084 


0.7219 


-0.1657 


0.3557 


2.4780 


2.00 


0.40521 


0.00186 


0.04959 


0.07074 


0.7169 


-0.1507 


0.3511 


2.4949 



to which the [L/M] Pade approximation is applied at 
the point 6 = 1. To perform the analytical continua- 
tion, the Pade approximant of [L/l] type may be used 
which is known to provide rather good results for various 
Landau- Wilson models (see, e.g., Refs. |H and0). The 
property of preserving the symmetry of a system during 
application of the Pade approximation by the 9 method, 
as in Ref. l43l . has become important for multivertices 
models. We used the [2/1] approximant to calculate the 
P functions in the two-loop approximation. 

However, analysis of the series coefficients for the /3 W 
function has shown that the summation of this series is 
fairly poor, which resulted in the absence of FP's with 
w* 7^ for a < 2.93. Dorogovtse\«2i found the symmetry 
of the scaling function for the WH model in relation to 
the transformation (u, v, w) — > (it, v, v + w), which gives 
the possibility of investigating the problem of the exis- 
tence of FP's with w* 7^ in the variables (it, v,v + w). 
In this case, our investigations carried out in RefQJ have 
shown the existence of FP's with w* ^ in the whole 
region where the parameter a changes. 

We have found two classes of FP's with g* = for 
an incompressible Ising model and with g* > for the 
compressible Ising model. If a possibility of realization of 
multicritical behavior in system for u* = u* — 2g* 2 = 
or for g* 2 = g* 2 — p,* = does not consider in this paper 
(for critical behavior 4 * 3 p,* = and g* = g*), then for 
either of the two classes there are three types of FP's in 
the physical region of parameter space, u*,g* 2 ,v*, and 
v* + w* > for different values of a. Type I corresponds 
to the FP of a pure system (u* ^ 0,v*,w* = 0), type II 
is a SR-disorder FP (u*,v* ^ 0, w* = 0), and type III 
corresponds to LR-disorder FP (u*,v*,w* ^ 0). 

The type of critical behavior of this disordered system 
for each value of a is determined by the stability of the 
corresponding FP. The stability properties of the FP's 
are controlled by the eigenvalues Ai of the matrix 



computed at the given FP: a FP is stable if all eigenvalues 
Xi are positive. If some of the eigenvalues A; are complex 



numbers, then the real parts of these eigenvalues must 
be positive for stable FP. 

Our calculations showed that FP's which belong to the 
class of incompressible Ising model with g* — are unsta- 
ble with respect to the influence of elastic deformations. 
For class of FP's corresponding to the compressible Ising 
model with g* ^ the values of the stable FP's ob- 
tained for 2 < a < 3 are presented in Table IIVI As 
one can see from this Table, for the compressible Ising 
model the LR-disorder FP is stable for values of a in the 
whole investigated range. The additional calculations for 
3 < a < 4 have shown that only SR-disorder FP is stable 
in this range. For a = 3 FP values for vertices u and 
g(k) are equal, u* = 0.26482 and v* + w* = 0.03448, 
and correspond to the SR-disorder FP for compressible 
Ising model, although w* ^ 0. Similarly for a = 3 the 
LR disorder is marginal, and the critical behavior of the 
compressible Ising model with LR-correlated disorder, as 
that of the SR-disordered compressible Ising model, is 
characterized by the same critical exponents (Table IIV[) . 

We have calculated the static critical exponents rj and 
v for the compressible Ising model with LR-correlated 
disorder (Table IIV[) , received from the resummed by the 
generalized Pade-Borel method 7 functions in the corre- 
sponding stable FP's: 

V = l<p{u* ,g\v* ,w*), (31) 

V = [2 + l!p 2{u\~g\v\w*)- 1 ^u\~g\v* 1 w*)]- 1 . 

Values of the specific-heat exponent a and the exponent 
P for magnetization which are required for calculation of 
ultrasound characteristics were determined from scaling 
relations a — 2 — dv and (3 — (d — 2 + 77)^/2 for 2 < 
a < 3 and presented in Table IIVI The values of the 
dynamic critical exponent z also presented in Table IIVI 
for 2 < a < 3 were taken from our paper^ where the 
critical dynamics of 3D disordered Ising model with LR- 
correlated disorder was considered for purely relaxational 
model A. It is caused by fact4£ that the coupling of the 
order parameter with elastic deformations is irrelevant for 
the relaxational critical properties of the order parameter 
in disordered compressible Ising model with a negative 
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FIG. 2: (a) Scaling functions for the critical sound attenua- 
tion 4>{y) an d (b) sound velocity dispersion f(y) in a double- 
logarithmic plot for disordered systems characterized by dif- 
ferent values of the correlation parameter a at T > T c and 
T < T c [4>o = 4>(0), fo = /(0)]. 



hydrodynamic region, where y ~ w£ z ~ (<z£) z "C 1, and a 
critical region y ~ w£ z 1 , which determines the behav- 
ior of the system near the phase transition temperature 
(r = (T — T c )/T c <C 1). We have seen from these curves 
that the correlation properties of a structural disorder 
does not affect the behavior of the scaling functions <j>(y) 
in the hydrodynamic region with y <C 1. However, LR- 
disorder begins to manifest itself in the crossover region 
1CP 1 < y < 10 and it has a drastic effect in the critical 
region with y > 10 (T -> T c ). 

As follows from Eqs. (fbl|) and (f2"Tj) , the attenuation 
coefficient can be expressed as 

o>,t)~ W 2 t" (32) 

and, using Eqs. (fT5)) and (|2^|) . we can write the relation 
for the sound velocity dispersion in the form 

cW)-c 2 (0,T)~T-"[/(y)-/(0)]. (33) 

The results of the calculation of the asymptotic de- 
pendences of the attenuation coefficient and the sound 
velocity dispersion for the critical and hydrodynamic re- 
gions are given in the Table [V] The characteristics of 
their frequency and temperature dependences were de- 
termined in the range 10~ 3 < y < 10" 1 for the hydro- 
dynamic regime and in the range 10 < y < 10 3 for the 
critical regime. Note that, according to,— the real tem- 
perature range 10~ 3 < r < 10 _1 in ultrasonic studies of 
phase transitions corresponds to the range 1 < y < 10 2 , 
i.e., it covers the crossover region and the beginning of 
the critical region (precritical regime). 

As it follows from the Table [V] the increase in cor- 
relation effects for disorder characterized by decrease in 
parameter a values leads to systematical increase in the 



specific-heat exponent a. 



IV. ANALYSIS OF RESULTS FOR 
ULTRASOUND CRITICAL CHARACTERISTICS 
AND CONCLUSIONS 

In this work, we employed the Pade-Borel approxima- 
tion of the summation of asymptotic series for calculation 
of the ultrasound scaling functions <j>(y) in (f2"4"| using ap- 
proximant [1/1]. The short series for scaling functions 
4>{y) were summed on the values of vertices u, v, and w 
in LR-disorder FP for 2 < a < 3 with variable y which 
was changed in the interval from 10~ 3 to 10 3 by steps 
with Ay = 0.1. 

The resultant behavior of the dynamic scaling func- 
tions (f>(y) and f(y) for individual values of a is shown in 
Figs.JHa) andf^b) on a log-log scale. Depending on the 
interval of changing variable y, the following asymptotic 
regions can be distinguished in the behavior of 4>(y): a 




-0.020 -0.015 -0.010 -0.005 0.000 0.005 0.010 0.015 



FIG. 3: Temperature dependence of the attenuation coeffi- 
cient a(uj, t) calculated for systems with LR-correlated dis- 
order characterized by different values of the correlation pa- 
rameter a at B — 0.3 and cu/Fo = 0.0015. 
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TABLE V: Asymptotic behavior of the sound attenuation coefficient a(c, r) in the critical, precritical, and hydrodynamic 
regimes for pure and disordered systems 



T <T C 


T > T c 


System 


Hydrodynamic 


Precritical 


Critical 


Hydro dy 


namic 


Precritical 


Critical 


Pure 




1.38 


W 1.08 T -0.21 


w 0.98 T -0.08 




1.38 


w 1.20 T -0.37 


W 1.05 T -0.17 


a = 3.0 




1.43 


w 1.21 T -0.24 


W 1.11 T -0.10 




1.43 


w 1.37 T -0.49 


w 1.21 T -0.24 


a = 2.9 


c 2 - 00 r- 


1.43 


w 1.30 T -0.27 


W 1.21 T -0.12 




1.43 


w 1.46 r -0.54 


w 1.29 T -0.28 


a = 2.8 




1.44 


w l-30 r -0.28 


W 1.21 T -0.13 




1.44 


w 1.46 r -0.55 


w 1.29 T -0.29 


a = 2.7 


c 2 - 00 r- 


1.47 


1 30 —0 29 
c r 


1 21 —0 13 


2 00 — 


1.47 


1 46 —0 56 


1 29 —0 29 
c r 


a = 2.6 


^•°°r- 


1.50 


w 1.30 T -0.30 


w 1.21 r -0.14 




1.50 


w 1.46 T -0.57 


w 1.30 T -0.30 


a = 2.5 


^•°°r" 


1.53 


w 1.30 T -0.32 


cl.21r-0.15 




1.53 


w 1.46 T -0.59 


w 1.30 T -0.30 


a = 2.4 


^•°°r- 


1.55 


w l-30 T -0.34 


W 1.21 T -0.16 




1.55 


w 1.47 T -0.63 


w l-30 T -0.32 


a = 2.3 




1.58 


w 1.32 T -0.37 


W 1.21 T -0.17 




1.58 


W 1.51 T -0.71 


w l-30 T -0.34 


a = 2.2 




1.60 


w l-32 T -0.39 


w 1.22 r -0.19 




1.60 


w 1.53 T -0.76 


w l-31 T -0.37 


a = 2.1 




1.62 


w 1.33 T -0.42 


W 1.23 T -0.21 




1.62 


w 1.56 T -0.83 


W 1.32 T -0.41 


a = 2.0 




1.64 


w l-34 T -0.45 


w 1.24 T -0.23 


J- 00 r- 


1.64 


w 1.60 T -0.92 


w l-33 T -0.44 



TABLE VI: Asymptotic behavior of the sound velocity dispersion c 2 (lu, t) — c 2 (0, r) in the critical, precritical, and hydrodynamic 
regimes for pure and disordered systems 



T < T c 


T>T C 


System 


Hydrodynamic 


Precritical 


Critical 


Hydro dynamic 


Precritical 


Critical 


Pure 


^2.00 r -2.65 


^0.30 r -0.49 


c°.Hr 


-0.25 




r -2.65 


W 1.08 T -1.48 


^0.34 r -0.54 


a = 3.0 


^2.00 r -2.93 


^0.38 r -0.51 


c°- 25 r 


-0.30 




r -2.93 


W 1.17 T -1.68 


^0.46 r -0.62 


a = 2.9 


^2.00 r -3.05 


^0.46 r -0.54 


c°- 32 r 


-0.33 


J™ 


r -3.05 


W 1.22 T -1.77 


w 0.57 r -0.72 


a = 2.8 


w 2.00 r -3.11 


w 0.46 T -0.55 


c°- 33 r 


-0.34 


a, 2 ' 00 


r -3.11 


W 1.32 T -1.98 


^0.61 r -0.80 


a = 2.7 


^2.00 r -3.16 


^0.46 r -0.56 


c°- 33 r 


-0.35 




r -3.16 


W 1.42 T -2.18 


w 0.65 T -0.88 


a = 2.6 


^2.00 r -3.21 


w 0.46 T -0.58 


c°- 33 r 


-0.36 


a, 2 ' 00 


r -3.21 


w 1.49 T -2.34 


^0.66 r -0.92 


a = 2.5 


^2.00 T -3.25 


w 0.46 T -0.60 


c°- 33 r 


-0.37 




r -3.25 


w 1.58 T -2.53 


w 0.66 r -0.94 


a = 2.4 


^2.00 r ~3.30 


^0.46 r -0.62 


^■ 33 r 


-0.38 


„ 2 - 00 


r -3.30 


w 1.62 T -2.64 


^0.67 r -0.97 


a = 2.3 


^2.00 T -3.34 


^0.47 r -0.65 


c°- 33 r 


-0.39 




r -3.34 


w 1.66 T -2.75 


W 0.67 T -1.00 


a = 2.2 


^2.00 r -3.38 


w 0.50 T -0.70 


c°- 34 r 


-0.39 


a; 2 ' 00 


r -3.38 


w 1.72 T -2.89 


W 0.68 T -1.03 


a = 2.1 


^2.00 r -3.41 


^0.51 r -0.74 


c°- 35 r 


-0.41 


a, 2 ' 00 


T -3.41 


w 1.76 T -2.98 


W 0.70 T -1.09 


a = 2.0 


w 2.00 r -3.43 


w 0.57 T -0.87 


c°- 37 r 


-0.51 


a; 2 - 00 


r -3.43 


w 1.80 T -3.06 


W 0.75 T -1.18 



values of exponents fci, and fcf* , which we can in- 
troduce for description of the critical anomalies of the 
frequency and temperature dependences of the atten- 
uation coefficient a ~ ' ' r^ k ^ ' and the sound ve- 
locity dispersion c 2 (lu,t) — c 2 (0,t) ~ uj k ^ r~ k ^ . So, 
for the attenuation coefficient the values of exponents 
ki a) = 1.21 and k ( T a) = 0.24 for T > T c and fc£ a) = 1.11 
and kr = 0.10 for T < T c in the case of point-like mi- 
correlated defects characterized by parameter a > 3.0, 
fc£ Q) = 1.29 and k ( T a) = 0.29 for T > T c and ki a) = 1.21 
and kr^ = 0.13 for T < T c in the case of complex struc- 
ture defects characterized, for example, by value of the 
parameter a = 2.7, and k^ = 1.30 and k^ = 0.32 for 
T > T c and ki a) = 1.21 and = 0.16 for T < T c in 
the case of extended linear defects characterized by value 



of the parameter a = 2.4. The comparison of these expo- 
nent values with the values of fcS = 1.05 and fef*^ = 0.17 
for T > T c and fc£ a) = 0.98 and k ( T a) = 0.08 for T < T c 
calculated in our paper^i for a pure Ising-like systems 
shows the strong influence of disorder and its correlation 
effects on the frequency and temperature dependences of 
the attenuation coefficient in the vicinity of the critical 
point. 

We must note that the increase in the attenuation co- 
efficient with increase in LR correlations for disorder as 
the critical temperature is approached is expected to be 
stronger than for systems with uncorrelated disorder or 
for the pure systems even in the hydrodynamic region. 
At the same time in the critical region, the systems with 
LR-correlated disorder should exhibit stronger both fre- 
quency and temperature dependences of the attenuation 
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coefficient than those for systems with uncorrelated dis- 
order or for pure systems. 

The similar conclusions we can make in relation to the 
frequency and temperature dependences of the sound ve- 
locity dispersion in the critical range on basis of those 

(c) (c) 

values of exponents ku and fcr which are given in Ta- 
ble ED and the values of ki c) = 0.34 and k [ T c) = 0.54 
for T > T c and k^ c) = 0.11 and k { T c} = 0.25 for T < T c 
calculated in paper— for pure systems. 

These conclusions are supported by the model repre- 
sentation of the results of the numerical calculations of 
the attenuation coefficient (Fig. [3|) for systems with LR- 
correlated disorder for different values of the parameter a 
performed at B — 0.3 and uj/Tq = 0.0015. These values 
were determined in Ref. |2l| when we compared the cal- 
culated temperature dependence of the attenuation co- 
efficient and the results of experimental studies of pure 
FeF2 samples^ which demonstrate Ising-like behavior in 
the critical region. 

A particularly important result of our investigation 
consists in the predicted manifestation of the dynamical 



effects of structural defects on anomalous ultrasound at- 
tenuation and dispersion over a wider temperature range 
near the critical temperature (already in the hydrody- 
namic region) in comparison with other experimental 
methods^ which require a narrow temperature range (of 
about t ~ 10~ 4 ) to be studied for revealing these effects. 
Thus, the results obtained in this paper can serve as a ref- 
erence point for purposeful experimental investigations of 
the dynamical effects of structural defects and their cor- 
relation properties on the critical behavior of solids using 
acoustic methods via the detection of the peculiarities of 
structural defects influence on the frequency and tem- 
perature dependences of the ultrasound attenuation and 
dispersion. 
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